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ABSTRACT 

In  a  seminal  paper  [1]  published  in  2000,  two  algorithmic  ver¬ 
sions  of  the  multichannel  parametric  adaptive  matched  filter 
(PAMF)  applied  to  space-time  adaptive  processing  (STAP) 
in  an  airborne  radar  application  were  shown  to  achieve  su¬ 
perior  test  detection  statistics  over  the  conventional  adaptive 
matched  filter  (AMF),  which  uses  a  non-parametric  approach 
to  estimate  the  detection  weight  vector.  In  fact,  the  perfor¬ 
mance  of  the  PAMF  approach  is  very  close  to  the  ideal  matched 
filter  (MF)  detection  statistics  under  exactly  known  covari¬ 
ance  (the  clairvoyant  case).  Improved  versions  of  the  two 
original  multichannel  PAMF  algorithms  [5],  one  new  mul¬ 
tichannel  PAMF  algorithm  [6],  and  a  new  two-dimensional 
(2D)  PAMF  algorithm  [all  four  with  fast  computational  im¬ 
plementations]  have  been  recently  developed.  In  this  paper, 
we  provide  the  outline  of  the  new  2D  parametric  algorithm 
and  summarize  the  detection  performance  of  3  of  the  4  new 
PAMF  algorithms  with  actual  Multi-Channel  Airborne  Radar 
Measurement  (MCARM)  data.  In  all  cases,  the  performance 
is  at  least  comparable  to,  and  in  some  cases  superior  to,  the 
original  multichannel  PAMF  algorithms  presented  in  [1],  while 
also  achieving  computational  savings  over  the  originals. 

1.  INTRODUCTION 

In  STAP  [2],  the  output  y  =  wTx  of  a  linear  combiner  (space- 
time  filter)  operating  on  the  radar  system  sensor  array  is  ex¬ 
pressed  in  terms  of  an  inner  product  of  a  multi-channel  (MC) 
weight  vector  w  of  dimension  M  x  1  and  a  comparable  di¬ 
mension  MC  data  vector  from  the  array.  It  is  well  known  that 
the  optimal  (adaptive)  weights  that  minimize  the  variance  of 
the  output  y  while  passing  a  signal  from  a  preferred  steering 
vector  direction,  represented  by  the  vector  e,  is 

w  =  R  xe  (1) 
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in  which  the  MC  covariance  matrix  R,/  is  based  on  the  second- 
order  statistics  of  the  disturbance  (no  target),  which  consists 
of  terrain  clutter  noise  plus  interferences  plus  receiver  white 
noise.  A  matched  filter  (MF)  detection  statistic  [1]  is  evalu¬ 
ated  by  comparing  the  value  of 
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in  which  e  is  steering  vector  and  x  is  the  current  measurement 
vector,  to  a  threshold  that  sets  the  probability  of  false  alarm. 
In  practice,  the  adaptive  matched  filter  (AMF)  detection  sta¬ 
tistic 
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is  evaluated  by  taking  the  inverse  of  the  estimated  covariance 
matrix 


Rd  =  (4) 

r=  1 

over  R  measured  realizations  of  the  MC  data  vector  x  of  di¬ 
mension  M.  This  is  usually  very  computationally  intensive 
if  the  dimension  of  the  M- channel  system,  which  is  typically 
the  product  of  the  number  of  array  channels  with  the  num¬ 
ber  of  pulses  processed,  is  in  the  hundreds.  As  an  alterna¬ 
tive  to  the  AMF,  [1]  and  [5]  discovered  that  multi-channel 
and  two-dimensional  parametric  estimation  approaches  could 
(1)  reduce  the  computational  requirements  by  reducing  the 
covariance  matrix  dimension  while  (2)  improving  the  detec¬ 
tion  performance  over  the  traditional  AMF  approach.  Both 
enhancements  were  achieved  by  replacing  the  inverse  of  the 
estimated  covariance  matrix  of  the  disturbance  with  the  order 
p  «  M  parametric  alternative 
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in  which  the  block  triangular  Toeplitz  matrix  structures,  de¬ 
veloped  in  [6],  are  functions  of  the  forward  MC  linear  predic- 
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tion  error  of  model  order  p  and  realization  r 

p 

e“[n;  r]  =  x[n;  r]  +  ^  Ap[k]x[n  -  fc;  r],  (6) 

/c=l 

the  backward  MC  linear  prediction  error  of  similar  model  or¬ 
der  p 

p 

ebp[n;  r\  =  x[n  - p;  r)  +  ^  Bp[fc]x[n  -p+k;r\,  (7) 

fc= l 

and  the  forward  and  backward  linear  prediction  error  vari¬ 
ances,  which  are  respectively 


2.  2D  LEAST-SQUARES  LINEAR  PREDICTION 
PARAMETER  ESTIMATION 

The  2-D  forward  linear  prediction  error  (LPE)  of  model  order 
p'2  for  the  2-D  signal  a;[nl,  n2;  r]  and  dimension  (pi  +  1)  x 
(X 1  —  pi)  is  defined  as 

P2 

ep2[™;  rl  =  rl  +  Ap2[fc]X[n  -  fc;  r]  (12) 

/c=l 

overp2+l  <  n2  <  N2,  in  which  AP2  [k\  are  the  2-D  forward 
linear  prediction  parameters  of  dimension  (pl  +  1)  x  (pl  +  1). 
This  may  be  expressed  in  vector  inner  product  form  as 


pp  =  £  {e?WKW)"}>pJ  =  £  {epM(epNr}  •  (8) 


The  MC  parameter  matrices  Ap[fe],  Bp[fc],  Pp,  and  Pp  are  es¬ 
timated  by  fast  computational  algorithms,  as  discussed  in  [3] 
and  [6],  One  algorithm  is  based  on  the  MC  harmonic-mean 
least-squares-based  lattice  algorithm  developed  by  Nuttall- 
Strand  [3]  [reference  as  MC-LLS]  and  the  other  is  based  on 
the  MC  geometric -mean  non-least-squares  lattice  algorithm 
developed  by  Viera-Morf  [3]  [reference  as  MC-LGM],  The 
resulting  parametric  adaptive  matched  filter  [PAMF]  detec¬ 
tion  statistic  for  these  two  algorithms  takes  the  form 


ApAMF-MF 
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(9) 


which  uses  the  parametric  inversion  formula  Eq.  5,  and  N  is 
the  number  of  pulses  processed.  An  exact  least-squares-based 
MC  covariance  case  of  MC  linear  prediction  uses  yet  another 
form  for  the  inverse  matrix,  specifically 
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in  which  the  various  block  triangular  Toeplitz  matrix  terms 
come  out  of  the  fast  computational  algorithm  discussed  in  [3], 
The  inverse  Eq.  10  is  used  in  lieu  of  If  ,/  in  the  PAMF  sta¬ 
tistic.  A  fourth  scheme,  covered  in  this  paper,  is  based  on 
two-dimensional  (2-D)  parametric  estimation  using  the  2-D 
least-squares-based  lattice  algorithm  [4].  The  specifics  of  the 
inverse  are  found  in  the  next  section  of  this  paper.  The  result¬ 
ing  2D  PAMF  detection  statistic  will  take  the  form 


ApAMF-2D  = 


E-=r9Et7e[j, 


[j> 


E7i  E„=1P e[j,n]«R.7e[j, 


(11) 


in  which  two  parametric  orders  p  [the  pulse  dimension]  and  q 
[the  array  element  dimension]  will  be  required. 


ep2[n;  r]  =  (I  Ap2)  Xp2[n;  r]  (13) 

for  which  Ap2  =  (AP2[1] . . .  AP2[p2])  is  a  1  x  p2  block  row 
vector  [scalar  dimension  (pi  +  1)  x  (pi  +  l)p2]  of  the  2-D 


forward  LP  parameters  and  X[n; 
/  x[pl,  k;  r] 

x[N  1  —  pi  —  1,  k;  r] 

V  x[0,k;r] 


rj  = 
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is  a  (pi  +  1)  x  (XI  —  pi)  data  matrix  and 

X[n;  r] 


2£p2[™;r]  = 


(15) 


X[n  —  p2;  r] 


is  a  (pi  +  l)(p2  +  1)  x  (XI  —  pi)  block  column  vector  of 
data  matrices.  Similarly,  the  2-D  backward  linear  prediction 
error  (LPE)  is  defined  as 


pi 

ep2[n;r]  =  X[n  —  p2;  r]  +  ^  Bp2[fc]X[n  —  p2  +  k\  r \ 

fc= l 

=  (Bp2  l)Xp2[n;r] 

for  which  Bp2  =  (Bp2[p2] . . .  Bp2[l])  is  a  1  x  p2  block  row 
vector  [scalar  dimension  (pi  +  1)  x  (pi  +  l)p2]  of  the  2-D 
backward  LP  parameters  (note  that  the  backward  LP  parame¬ 
ter  vector  has  reversed  indexing  relative  to  the  forward  LP 
parameter  vector).  The  (pi  +  1)  x  (pi  +  1)  2-D  forward  LPE 
variance  is  both  defined  and  expressed  as 

Pp2  =  £  {ep2(ep2)H}  =  (I  ap2)  rp2  (  j  (16) 

in  which  Rp2  is  the  block  (p2  +  1)  x  (p2  +  1)  2-D  double 
Toeplitz  autocovariance  matrix  [it  is  also  hermitian  symmet¬ 
ric]  with  a  scalar  dimension  of  (pi  +  l)(p2  +  1)  x  (pi  + 


(26) 


l)(p2  +  1)  with  structure 
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composed  of  the  stationary  (pl+1)  x  (pl+1)  scalar-dimensioned 
2-D  autocovariance  matrix  elements  R[fc]  = 
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B  A 

The  block  diagonal  matrix  Pp2  is  similar  to  Pp2,  except  the 
block  diagonal  elements  are  (JPp2J)_1. 


\r[-pl,k]  r[—(pl  —  1),  k]  •••  r[0,fc]  / 

(18) 

Note  that  Rp2  is  double  Toeplitz  (Toeplitz  at  the  block  level 
and  also  the  block  matrix  elements  themselves  are  Toeplitz); 
this  property  results  in  the  relationship  Rp2  =  JR*2J  and 
Rp2  =  Rp2.  In  a  similar  manner,  the  2-D  (pi  + 1)  x  (pi  + 1) 
backward  LPE  variance  is  given  by 

Pp2  =  £  {ebp2(ebp2)H}  =  (Bp2  I)  Rp2  (  -f  )  •  (19) 

Assuming  the  initialization  Po  =  R[0],  the  four  steps  of  the 
2-D  Levinson  recursion,  that  solves  for  the  LP  parameters  Ak 
and  2-D  LPE  variance  P*.  over  k  =  1, ...  ,p,  are 


/  R[fc]  \ 

Ak  =  (lAk_1)  :  =JAP 

(20) 

V  m ) 

rfc  = -Afc  (JP^J)"1 

(21) 

(I  A*)  =  (I  Afe-1  0)  +  Tk  (0  JA^J  I) 

(22) 

pfe  =  (i-rfc(jr*j))pfe_i. 

(23) 

This  will  reduce  to  1-D  Levinson  recursion  when  pi  =  0; 


3.  PARAMETRIC  INVERSE  OF  2D  DOUBLY 
TOEPLITZ  COVARIANCE  MATRIX 


The  form  of  the  doubly  Toeplitz  block  inverse  [4]  which  is 
useful  to  this  project,  and  which  also  depends  only  on  the  last 
parameter  order  computed,  is 


R— 1  _  A  H  P  t  A  _  TdH  -p-B  [_> 

p2  —  — p2— p2±Lp2  J=*p2£-P2l=*p2 


(24) 


in  which  the  block  triangular  Toeplitz  matrices  are  given  as 
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4.  PERFORMANCE  ANALYSIS  USING  ACTUAL 
MCARM  RADAR  DATA 

For  details  on  the  motivation  and  structure  of  the  MCARM 
radar  data,  the  reader  is  referred  to  Section  C  of  reference  [1], 
A  key  issue  in  performance  testing,  specifically  the  probabil¬ 
ity  of  detection  Pd,  is  the  amount  of  secondary  (training)  data 
[small  vs  large  sample  sizes]  required  to  achieve  a  reasonable 
level  of  detection  performance  in  a  constant  false  alarm  rate 
Pfa  (CFAR)  scenario.  Three  scenarios  that  inject  artificial 
target  signals  into  the  actual  MCARM  data  were  considered, 
leading  to  eight  plot  sets  of  test  statistic  magnitude  vs  range 
bin  index  for  various  multi-channel  and  2D  algorithm  order 
selections.  The  MCARM  database  used  for  this  analysis  ex¬ 
tracted  (from  acquisition  575,  flight  5)  one  elevation  channel, 
4  azimuth  channels,  32  pulses,  and  inclusive  range  bins  (RBs) 
142-469  from  each  pulse.  In  test  case  1,  a  fixed  window  of 
training  data  utilized  256  RBs  (142-269  and  341-468),  a  pair 
of  simulated  primary  targets  with  an  input  SINR  of  -30  dB 
were  inserted  into  RBs  291  and  293,  and  the  appropriate  de¬ 
tection  test  statistic  (either  MC  Eq.  9  or  2D  Eq.  1 1)  generated 
from  RBs  270-340  for  plotting.  In  test  case  2,  a  moving  win¬ 
dow  of  training  data  consisted  of  11  adjacent  RBs  centered 
sequentially  over  RBs  270-340,  with  the  same  simulated  pri¬ 
mary  targets  inserted  into  RBs  291  and  293,  and  the  test  sta¬ 
tistics  then  generated.  The  three  center  RBs  in  the  moving 
window  were  not  used  in  the  estimate  of  the  inverse  covari¬ 
ance,  leaving  only  a  short  8  RB  secondary  data  sequence  (in 
contrast  to  the  long  256  RB  test  case  1).  The  test  case  3  is 
identical  to  test  case  1,  except  four  secondary  disturbance  tar¬ 
gets  were  artifically  inserted  into  RBs  238,  269,  373,  and  400, 
each  with  an  input  SINR  of  -10  dB  and  the  same  steering  vec¬ 
tors  as  the  primary  targets,  in  order  to  evaluate  the  effects  of 
targets  in  the  secondary  data. 

Figure  1  contains  extracted  test  statistic  plots  from  [1]  at  a 
fixed  parametric  order  that  demonstrate  the  MC  lattice-based 
algorithm  performance  (labeled  PAMF-SN)  and  the  MC  co- 
variance  least-squares-based  algorithm  performance  (labeled 
PAMF-LS)  for  the  three  test  case  conditions.  In  all  cases,  the 
reference  found  that  both  parametric  MC  algorithms  yielded 


TEST  CASE  1 


TEST  CASE  2 


TEST  CASE  3 


Fig.  1.  PAMF-SN  and  PAMF-LS  test  statistics  from  [1]  for  3 
test  cases  to  compare  to  3  new  MC  and  2D  algorithms. 


Fig.  2.  Test  Case  1  (fixed-window  filter  and  two  primary  tar¬ 
gets):  PAMF-SN  test  statistics  of  multi-channel  lattice  algo¬ 
rithm  for  MC  orders  1  to  6. 


better  detection  statistics  than  the  conventional  non-parametric 
adaptive  matched  filter  (AMF)  algorithm  approach.  Figures 
2-8  are  test  statistic  plots  of  the  two  revised  MC  SN  and  LS 
algorithms  and  the  new  2D  algorithm  over  a  range  of  paramet¬ 
ric  orders  from  1  to  6.  The  reader  will  notice  that  the  perfor¬ 
mance  is  not  very  sensitive  to  model  order,  so  comparison  will 
be  made  with  Figure  1  performance  at  the  order  indicated  in 
that  figure.  All  plots  were  normalized  to  set  the  mean  value  to 
0  dB  for  common  referencing.  Of  interest  as  indicators  of  the 
quality  of  the  test  statistics  are  two  number  pairs:  (1)  the  max 
target  peak  above  the  0  dB  reference  level  and  (2)  the  differ¬ 
ence  between  the  target  peak  value  and  the  highest  non-target 
peak  value  (desire  large  values  for  both  of  these  numbers).  In 
Figures  2  and  3  at  order  3  for  test  case  1,  the  new  SN  and 
LS  MC  algorithms  produced  number  pairs  of  (23dB,13dB) 
and  (22dB,14dB)  respectively,  compared  to  (25dB,16dB)  and 


Fig.  3.  Test  Case  1  (large  fixed-window  filter  and  two  primary 
targets):  PAMF-LS  test  statistics  of  multi-channel  covariance 
least  squares  algorithm  for  MC  orders  1  to  6. 


Fig.  4.  Test  Case  2  (short  moving-window  filter  and  two  pri¬ 
mary  targets):  PAMF-SN  test  statistics  of  multi-channel  lat¬ 
tice  algorithm  for  MC  orders  1  to  6. 


(31dB,19dB)  in  Figure  l(a)(b).  In  Figures  4  and  5  at  order  3 
for  test  case  2,  the  improved  SN  and  LS  MC  algorithms  pro¬ 
duced  number  pairs  of  (32dB,22dB)  and  (23dB,13dB)  respec¬ 
tively,  compared  to  (25dB,16dB)  and  (22dB,lldB)  in  Figure 
l(c)(d).  In  Figures  6,  7,  and  8  for  test  case  3,  the  new  SN  and 
LS  MC  algorithms  and  the  new  2D  algorithm  produced  num¬ 
ber  pairs  of  (29dB,17dB),  (24dB,13dB),  and  (22dB,10dB)  re¬ 
spectively,  compared  to  (15dB,6dB)  and  (21dB,l  ldB)  for  the 
MC  cases  in  Figure  l(e)(f). 


Fig.  5.  Test  Case  2  (short  moving-window  filter  and  two  pri¬ 
mary  targets):  PAMF-LS  test  statistics  of  multi-channel  co- 
variance  least  squres  algorithm  for  MC  orders  1  to  6. 


Fig.  7.  Test  Case  3  (large  fixed-window  filter,  two  primary 
targets,  four  secondary  data  targets):  PAMF-LS  test  statistics 
of  multi-channel  lattice  algorithm  for  MC  orders  1  to  6. 


Fig.  6.  Test  Case  3  (large  fixed-window  filter,  two  primary 
targets,  four  secondary  data  targets):  PAMF-SN  test  statistics 
of  multi-channel  lattice  algorithm  for  MC  orders  1  to  6. 
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AIRBORNE  RADAR 
SURVEILLANCE  APPLICATION 


Uniform  Planar  Antenna  Array 


MCARM:  Multi-Channel  Airborne 

Radar  Measurement  Platform 


ULA:  Uniform  Line  Array 
UPA:  Uniform  Planar  Array 

Data  is  four-dimensional  (  4-D  ) 
J-22 
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DATA  ORGANIZATION  POSSIBILITIES 
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c)  •  Multi-Channel  (M-C)  +  “Realizations”  [  window  of  range  bins  (RBs)  ] 


x 


jn 


[k] 


Note:  Numbers  shown  are  for  MCARM  radar  data.  — 


A 


Conventional  multi-channel  STAP 
configuration  for  arbitrary  array 
geometry  or  non-uniform  phase 
centers  of  planar  array  elements 
[“channel”  is  a  paired  antenna 
element  position  and  a  pulse  number  ] 


to 

T» 

00 

«N 

II 


v 


-  k  =  1 


k =  K<360 


e)  •  Two-Dimensional  (2-D)  +  “Realizations”  (RBs) 

x\jh,ri\  k] 


Appropriate  for  uniform  line  arrays 
(  ULAs  )  [  or  a  single  row  of  an  UPA  ] 


Three-Dimensional  (3-D)  +  “Realizations”  (RBs) 

x[jhJv,n;  k] 


k  =  K 


k  =  1 


T-A 


Appropriate  for  uniform  planar 
arrays  (  UPAs  )  [  both  rows  &  cols  ] 


//////// 

k  =  1 


/// 
/// 
/// 


M  k  =  K 


Four-Dimensional  (4-D) 

x[jh,3v,n,k] 


All  dimensions  of  data  are  used;  no 
dimension  is  available  for  ensemble 
“realizations” 


ASAP2006-3 


MATCHED  FILTER  AND  DETECTION  STATISTIC 
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•  M-C  Linear  Combiner  /  Waveformer: 


y  =  W/7X 


H  means  complex  conjugate  transpose 
T  means  transpose 


uses  JN  (#  antenna  array  elements  x  #  pulses)  M-C  data  vector: 

—  r  Y 

X.  I  X  |  . . .  X  i  X  2  •  •  •  X  ^  X  jy  ■  ■  •  X  jy  I 

•  M-C  Waveformer  MSE  Weights  To  Create  Matched  Filter  w^p  =  R~^e 

in  which  Rr/  is  the  JN  x  JN  M-C  target-free  “disturbance”  covariance  matrix 
and  e  is  a  JNx'\  desired  steering  vector  of  the  receiver  array  (a  function 
of  wavelength,  spatial  spacings,  angles  of  arrival,  PRI,  target  velocity) 

•  M-C  Matched  Filter  Waveformer  Output  Usually  evaluated  for  values 

1  of  e  (normalized  Doppler  freqs. 
VMF  =  wMFx  =  e  R  ,  X  or  normalized  AOA  directions); 

a  M  =  256  to1024  is  nominal 

•  M-C  DETECTION  STATISTIC:  Normalized  Magnitude  of  Matched  Filter  Output 

eHK^1x2  i - ii - m  _ 

AMF  =  (e*R-e)  (x«R-x)  UU  =  ' 

•  DETECTION  STATISTIC  COMPARED  TO  A  Pr {fa}  THRESHOLD  FOR  DETECTION 
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ALGORITHM  DEVELOPMENT  OBJECTIVES 
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•  COVARIANCE  MATRIX  INVERSE  ESTIMATION  STRATEGIES  THAT  : 

•  Significantly  reduce  computational  requirements 

•  Inversion  of  covariance  matrix 

•  Evaluation  of  matched  filter  vector-matrix  products 

•  Reduce  secondary  (training;  RBs)  data  requirements 

•  Increase  detectability  for  same  or  less  data 


•  APPROACHES  THAT  USE  : 

•  Parametric  vs  non-parametric  estimation  techniques 

•  Least  square  error  (LSE)  vs  mean  square  error  (MSE) 

•  Primarily  multi-channel  (M-C)  structures;  also  try  2-D  structures 
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STRUCTURAL  FORMS  OF  COVARIANCE  MATRIX 
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•  Assume  uniform  PRI  (  equi-spaced  N )  and  ULA  (  J=Jh  and  V=1  ) 

•  Multi-channel  (M-C)  form  (uncalibrated  array) 


Block  Nx  A/Toeplitz  and  Hermitian  symmetry 
Element  J  x  J  dimension  without  special  structure 


N  blocks 


•  Two-dimensional  (2-D)  form  (calibrated  array) 


Block  Nx  A/Toeplitz  and  Hermitian  symmetry 


Element  J  x  J  Toeplitz  without  special  symmetry 


N  blocks 


ASAP2006-6 


CONVENTIONAL  NON-PARAMETRIC  APPROACH 
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•  Formulation  of  non-parametric  test  statistic 

•  Often  referenced  as  the  Adaptive  Matched  Filter  (AMF)  test  statistic 

•  Generate  covariance  matrix  estimate  from  range-indexed  data  vectors 

4=^1 ]x[k]x[kf  order  K(JN)2 

K  k= 1 

•  Compute  covariance  matrix  inversion  R~l  order  (  N) 

•  Form  MF  quadratic  product  terms  eH  Rdx ,  eH Rdxe ,  xhRdx 

for  a  range  of  M  steering  vector  e  values  using  x  data  vector  not  in  Rdl 

order  M  (JN)2 
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ALTERNATIVE  APPROACH  #1  OVERVIEW  : 
MSE-BASED  M-C  AR  PARAMETRIC  INVERSE 


OSUI 

Oregon  State 


UNIVERSITY 


•  Range-indexed  data  ->  M-C  autoregressive  (AR)  block  JxJ 
parameter  estimates  (several  M-C  algorithms) 


•  Express  inverse  in  terms  of  products  of  block  triangular  matrices 
formed  from  the  block  AR  parameters 


block  diagonal  matrix  with  inverses  of  error 
covariances  along  diagonal 


block  triangular  matrices  with  block  AR  parameters 


•  Quadratic  product  terms  exploit  block  triangular  matrices  to  reduce 
computational  evaluation.  Example: 


special  triangular  matrix-vector  products 
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ALTERNATIVE  APPROACH  #2  OVERVIEW  : 
LSE-BASED  M-C  LP  PARAMETRIC  INVERSE 
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Not  forced  to  use  block-Toeplitz  form  of  #1  MSE-based  approach 


•  Inverse  R/  replaced  by  inverse  of  block-data  matrix  product  Xd) 

X d  —  rectangular  block-Toeplitz  structure  formed  by  stacking  x[£]  data  vectors 

•  Range-indexed  data  ->  M-C  linear  prediction  (LP)  block  JxJ 
parameter  estimates  (two  M-C  LP  algorithms*) 


-i 


Express  inverse  in  terms  of  products  of  block  triangular  matrices 
formed  from  the  block  M-C  LP  parameters 


block  diagonal  matrix  with  inverses  of  least  squared 
error  sums  along  diagonal 

block  triangular  matrices  with  block  LP  parameters 


•  Quadratic  product  terms  exploit  block  triangular  matrices  to  reduce 
computational  evaluation.  Example: 


special  triangular  matrix-vector  products 


*  Marple:  Digital  Spectral  Analysis,  Prentice  Hall,  1987  (chap. 15)  and  two-volume  Signal  Analysis,  2008 
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1-D  MOTIVATIONAL  EXAMPLE  : 

MINIMUM  VARIANCE  SPECTRAL  ESTIMATION 
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Order  p  minimum  variance  (Capon)  spectral  estimators 


MSE-based  form  (Toeplitz  covariance  matrix) 

=  7 - -  for  — —  </<—*—  ,  T  =  sample  interval 


MV-MSE  ' 


IT 


IT 


f 


eH  (/)  =  [!  exp(-y2^/TJ  •••  exp(-j2;rpTj\  ,  R  = 


r[  0]  r*[  1] 
r[  1]  r[  0] 


r*[p] 


:  •.  •.  r[l] 

\r[p]  •••  r[!]  r[°]  J 

•  LSE-based  form  (non-Toeplitz  covariance-like  data  matrix  product) 
for  data  record  of  N  samples 

f x\p  + 1]  •••  x[N-p]  •••  x[N]  ^ 

(/)  =  - 1 - T -  •  X*  = 


R 


MV-LSE ' 


eH (f)(XH xy' e(f) 


V 


x[l] 


x[p  +  1] 


x[7V-/>] 


J 


•  Performance  test  case:  3  complex  sinusoids  in  additive  white  noise  (5  dB  SNR) 
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1-D  MOTIVATIONAL  EXAMPLE  : 
PERFORMANCE  COMPARISON 
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ESTIMATED  MSE-BASED  MINIMUM  VARIANCE  SPECTRUM 


Frequency  (Hertz) 


ESTIMATED  LSE-BASED  MINIMUM  VARIANCE  SPECTRUM 


-1000  -800  -600  -400  -200  0  200  400  600  800  1000 

Frequency  (Hertz) 


ESTIMATED  LSE-BASED  MINIMUM  VARIANCE  SPECTRUM 


-1000  -800  -600  -400  -200  0  200  400  600  800  1000 

Frequency  (Hertz) 


MSE-Based 
(Data  Correlation) 
Minimum  Variance 
Spectral  Estimate 

N=200  samples 


LSE-Based 
(LP  Correlation) 
Minimum  Variance 
Spectral  Estimate 

N=200  samples 


LSE-Based 
Minimum  Variance 
Spectral  Estimate 

N-14  samples 


All  spectral  estimates  above  are  order  p  =  6 
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1-D  MOTIVATIONAL  EXAMPLE: 

BASIS  FOR  MSE-BASED  FAST  ALGORITHM 


*  Covariance  matrix  R  isToeplitz 
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Non-Toeplitz  triangular  decomposition  (Cholesky) 
_  H 


f  Pt 
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1  a2[  1]  a^[  2] 
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0  0  1 


Note  that  all  orders  from 
1  to  p  are  used 


Toeplitz  triangular  decomposition  (Gohberg-Semencul) 
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Note  that  only  order 
p  is  used 


Obtain  AR  parameters  by  solution  of  Yule-Walker  equations  (Levinson) 
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1-D  MOTIVATIONAL  EXAMPLE: 
MSE-BASED  FAST  ALGORITHM 
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Substituting  parametric  inverse  covariance  matrix  yields 

1  1 

^  FFT  operation  ! 


^MV-MSE  CO 


eH  (f)Rp'e(f)  Yk-_.  ’PI*]  vtv(~j2njkTs ) 


t=~P 


Psi  coefficients  are  autocorrelations  of  AR  parameters  (  2  FFTs  ! ) 


¥[£]  =  ¥*[—£]  = 


i  p- 1 

~lL{p+l-k  -2 i)a  \k  +  i]a'[i]  0 <k<p 

Pp  i= 0 


0  k  >  p 

Many  AR  parameter  estimation  algorithms  are  available  : 


•  Levinson:  x  ->  r  ->  a 

•  Burg:  x  ->  a 

•  Fast  computational  algorithm  order:  Np  +  3p\og2  p 
Compare  to  direct  covariance  matrix  inversion:  J\J  2 


Musicus,  IEEE  Trans.  ASSP,  pp.  1333-1335,  October  1986 
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1-D  MOTIVATIONAL  EXAMPLE: 
LSE-BASED  FAST  ALGORITHM 


Fast  computational  structure  of  LSE-based  MV  spectral  estimator 

1  1 


R 


MV-LSE 


(/)  = 


eH  (/)  ( X”Xp  )_1  e(f)  VW  e*p(-j2xJkTs ) 


Based  on  following  4-term  triangular  Toeplitz  inverse  structure 

H _ _ H _ H _ _ H 


forward  LP  backward  LP  forward  time 

parameters  &  parameters  &  edge  LP  gain 

1st  sqs  error  sum  1st  sqs  error  sum  vector 


backward  time 
edge  LP  gain 
vector 


•  All  4  parameter  sets  are  obtained  from  the  fast  solution  algorithms*  of 
forward  and  backward  least  squares  linear  prediction 


•  Psi  parameter  'i'  is  the  sum  of  4  autocorrelations  of  the  4  parameters 
(compared  to  1  autocorrelation  in  the  prior  MSE  case,  where  b  —  a  ) 


•  Overall  computational  complexity  of  order  Np2  +  9p\og2p 


*  Marple:  Digital  Spectral  Analysis,  Prentice  Hall,  1987  (chap. 15)  and  two-volume  Signal  Analysis,  2008 
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1 -D  MOTIVATIONAL  EXAMPLE: 
CONCLUSIONS 
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•  LSE-based  minimum  variance  spectral  estimation 

•  Fast  computational  algorithm(s)  are  available  that  are 
competitive  with  MSE-based  versions 

•  Better  performance  is  achieved  over  MSE-based  versions 
(with  same  data  record  size) 

•  Similar  performance  is  achieved  as  that  of  MSE-based 
versions  (with  significantly  decreased  data  record  size) 
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BACKGROUND:  M-C  PARAMETRIC  MATCHED 
FILTER  AIRBORNE  RADAR  APPLICATION 
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•  Seminal  April  2000  IEEE  AES  Transaction  Paper* 

•  Replaced  non-parametric  estimate  of  covariance  inverse  with  parametric  est. 

•  Applied  non-Toeplitz  triangular  inversion  formula,  so  no  FFT  eval  of  MF  terms 

•  Used  M-C  extension  of  Burg  algorithm  by  Nuttall-Strand**  to  est.  M-C  AR  terms 

•  Statistical  averaging  over  range  of  disturbance-associated  data  range  vectors 
was  applied  outside  M-C  AR  parameter  estimation  (post-estimation) 

•  Summer  2006  AFOSR  Fellow  Study 

•  Applied  Toeplitz  triangular  inversion  formula,  so  fast  FFT  eval  of  MF  terms 

•  Modified  M-C  Nuttall-Strand  algorithm  to  average  range  inside  algorithm 

•  Applied  M-C  extension  of  Marple  fast  M-C  LP  algorithm***  to  estimate  M-C 
LSE-based  LP  parameters  to  use  in  the  inversion  formula 

•  Features  For  Comparison:  fast  computations  and  increased  performance 

*  I^^S^^<angaswamy,  Davis, ^hang,  Himed',  Michels,  vol.^ib,  pp7fc77-69iT^000 

**  Marple,  Digital  Spectral  Analysis  (1987)  and  ***  Signal  Analysis  (2008) 
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KEY  COMPUTATIONAL  COMPONENTS  OF  FAST 
M-C  MSE-BASED  PARAMETRIC  ALGORITHM 


Fast  M-C  MSE-based  AR  parameter  estimation  algorithm  follows 
similar  approach  as  1-D  motivational  example,  except  block  matrix 
structure  replaces  scalar  parameters 

H  H 

t 


dim  JN-^Jp 
[eg,  A/=256,  p=6] 


\ 


block  order  p2 


o p2Ar 


f i  4W 

A2[  2f! 

f 

0 

i 

4[i] 

vO 

0 

i  J 

V 

0 

0 

B\- 1 


(P2) 


o  b;  [2] 
o  o 

o  o 


b;i 2] 
0 


J 


Use  enhanced  Nuttall-Strand-Marple  M-C  AR  estimation  algorithm  to  get  both  forward 
and  backward  block  AR  parameters  (conjugate  relationship  does  not  hold  in  M-C) 


M-C  matched  filter  (MF)  evaluation  for  detection  statistic  exploits  inverse  structure 


three  block  FFTs, 
block  order  3 p  log2  p 


special  block  linear  convolution  products  (only  1  of  2  terms  shown) 
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EXPERIMENTAL  BASIS  FOR 
PARAMETRIC  ORDER 
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ACTUAL  MCARM  RADAR  DATA 


a> 

o 

CL 

TO 

+-» 

O 

CD 

Q. 

CO 


Multi-Channel  AR  Autospectrum  of  Array  Element  1  at  Order  6 


-  Doppler  0  Doppler  + 

Multi-Channel  AR  Autospectrum  of  Array  Element  3  at  Order  6 


Multi-Channel  AR  Autospectrum  of  Array  Element  2  at  Order  6 


Multi-Channel  AR  Autospectrum  of  Array  Element  4  at  Order  6 


SIMULATED  RADAR  DATA 

Multi-Channel  AR  Autospectmm  of  Array  Element  1  for  KJcdfeOiannel  AR  Autospectrum  of  Array  Element  2  for  Order  3 


Multi-Channel  AR  Autospectrum  of  Array  Element  3  for  BatteOiannel  AR  Autospectrum  of  Array  Element  4  for  Order  3 


Conditions: 

1.  Sensor  array  channels  1  to  4  shown  (out  of  22) 

2.  Order  6  MV  spectral  estimates  shown 
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KEY  COMPUTATIONAL  COMPONENTS  OF  FAST 
M-C  LSE-BASED  PARAMETRIC  ALGORITHM 
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•  Similar  to  M-C  MSE-based  block  parametric  algorithm,  except  have 
twice  as  many  block  triangular  components,  but  can  apply  block 
FFT  operations 


•  Based  on  following  4-term  block  triangular  Toeplitz  inverse  structure* 


forward  block 
LP  parameters  & 
1st  sqs  error  sum 


backward  block 
LP  parameters  & 
1st  sqs  error  sum 


forward  time 
edge  LP  block 
gain  vector 


backward  time 
edge  LP  block 
gain  vector 


dim  N->Jp 
[eg,  N= 256,  p= 6] 

block  order  p2 


special  block  linear  convolution  products  (only  1  of  4  terms  shown) 


*  Marple,  Digital  Spectral  Analysis  (1987)  and  Signal  Analysis  (2008) 


ASAP2006-19 


KEY  COMPUTATIONAL  COMPONENTS  OF  FAST 
2-D  MSE-BASED  PARAMETRIC  ALGORITHM 
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*  Fast  2-D  MSE-based  AR  parameter  estimation  algorithm  follows 
similar  approach  as  1-D  motivational  example,  except  block  matrix 
doubly  Toeplitz  structure  replaces  scalar  Toeplitz  parameters 


f(p2  r1 

0 

0 


0 


(P2) 

0 


dim  JN->PjPn 
[  eg,  J= 22,  A/=256, 
Pj=4,  PN= 6  ] 


block  order  pN2 , 
scalar  order  Pj2pN2 


•  Use  2-D  AR  estimation  algorithm  published  by  Marple  [IEEE  Signal  Processing 
Letters,  2000]  to  get  2-D  and  backward  block  AR  parameters 


2-D  matched  filter  (MF)  evaluation  for  detection  statistic  exploits  inverse  structure 


special  block  doubly  Toeplitz  linear  convolution  products  (only  1  of  2  terms  shown) 
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PERFORMANCE  COMPARISONS 
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•  SIMULATED  RADAR  ARRAY  DATA 

•  Motivated  by  MCARM  data  collection 

•  Six  measurements  conditions  investigated  (3  shown  here) 

•  Plot  probability  of  detection  vs  output  SINR  for  fixed  Pr{fa} 


•  ACTUAL  RADAR  ARRAY  DATA  (MCARM  public  data) 

•  Subset  of  total  available  data 

•  Three  test  cases  with  artificially  inserted  target  and  interference  signals 

•  Plot  MC  or  2D  detection  statistic  vs  range  (about  20  Km) 
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SIMULATION  PARAMETERS 
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•  Based  on  Ward’s  Clutter  Model 


-  Side-looking  GMTI  Scenario;  modeled  after  April  2000  IEEE  Paper 

z 

▲ 


•  MCARM-like  scenario 

•  40  dB  Clutter-to-Noise 
(CNR)  Ratio 

•  2  Jammers  (some  cases) 

•  20°  Crab  (some  cases) 

•  J=4  Elements 

•  N=32  Pulses 

•  Sample  support  (either): 

-  K=256  range  bins 

-  K=8  range  bins 
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SUMMARY  OF  SIMULATED  TEST  CONDITIONS 
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SIMULATED  TEST  CONDITIONS  CONSIDERED 


Large  training  set,  no  crabbing,  no  shared-angle  interferers 
Small  training  set,  no  crabbing,  no  shared-angle  interferes 
Large  training  set,  crabbing,  no  shared-angle  interferers 
Small  training  set,  crabbing,  no  shared-angle  interferers 
Large  training  set,  no  crabbing,  two  interferers 


Small  training  set,  no  crabbing,  two  interferers 


•  PERFORMANCE  MEASURE:  dB  offset  to  right  of  ideal  MF  clairvoyant  response  (shown 
on  next  slide)  at  Pd  of  0.5  (represents  roughly  the  performance  degradation  of  the  curve) 


•  PERFORMANCE  SUMMARY  TABLE 


M-C  Lst  Sqs  Lattice 

M-C  Geometr.  Lattice 

M-C  Lst  Sqs  Linear  Predict. 

2-D  Lst  Sqs  Lattice 

2006  Study 

2000  Study 

2006  Study 

2000  Study 

2006  Study 

2000  Study 

2006  Study 

2000  Study 

Test  Case  1 

-0.5  dB 

-1  dB 

-0.5  dB 

_  _ 

-0.5  dB 

-0.5  dB 

-0.5  dB 

mm 

Test  Case  2 

-7.5  dB 

-8  dB 

-0.7  dB 

— 

-0.5  dB 

-0.5  dB 

-1.2  dB 

— 

Test  Case  3 

-0.8  dB 

-1  dB 

-1.0  dB 

— 

-0.8  dB 

-0.8  dB 

-3.0  dB 

— 

Test  Case  4 

-3.5  dB 

-3.5  dB 

-1.0  dB 

— 

-0.8  dB 

-0.8  dB 

-3.0  dB 

— 

Test  Case  5 

-0.6  dB 

-0.6  dB 

-0.6  dB 

— 

-0.6  dB 

-0.6  dB 

-1.2  dB 

— 

Test  Case  6 

-2.5  dB 

-2.5  dB 

-1.2  dB 

- 

-0.6  dB 

-0.6  dB 

-1.5  dB 

- 

•  LstSqs  LP  (MC-LLS)  performs  most  robustly  with  little  variations  across  all  test  conditions 

•  Few  tenths  dB  improvement  of  improved  LstSqs  Lattice  algorithm  (MC-CLS  2006)  over 
MC-CLS  2000  algorithm  version  for  several  test  conditions 
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SIMULATED  TEST  CONDITION  #1 
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J= 4  Array  Elements,  N=32  Pulses,  K= 256  Fixed  Range  Bin  Window 
No  Crab  Angle,  No  Interferes  Added,  Order  o=3 

2006  Study  2000  Study 
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SIMULATED  TEST  CONDITION  #2 


OSUI 

Oregon  State 


UNIVERSITY 


J= 4  Array  Elements,  N=32  Pulses,  K=  8  Moving  Range  Bin  Window 
No  Crab  Angle,  No  Interferes  Added,  Order  p= 3 

2006  Study  2000  Study 
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SIMULATED  TEST  CONDITION  #6 
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J= 4  Array  Elements,  N=32  Pulses,  K=  8  Moving  Range  Bin  Window 
No  Crab  Angle,  2  Jammers  Added,  Order  p= 3 

2006  Study  2000  Study 
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CONCLUSIONS 
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•  Parametric  covariance  matrix  fast  inversion  algorithms  for  : 

•  Enhanced  M-C  MSE-based  (force  block  Toeplitz)  inversion  developed 

•  M-C  LSE-based  (non-block-Toeplitz)  inversion  introduced 

•  2-D  MSE-based  (force  doubly  Toeplitz)  inversion  developed 


•  Using  the  non-parametric  inversion  approach  as  the  baseline, 
the  relative  computational  requirements  and  performance  as 
averaged  over  the  3  simulations  shown  were  : 

•  M-C  MSE-based  algorithm  was  .7  x  baseline  in  computations  and 
-.5  dB  from  exact  MF,  vs  -2.8  dB  for  non-parametric  baseline 

•  M-C  LSE-based  algorithm  was  1.1  x  baseline  in  computations  and 
-.5  dB  from  exact  MF,  vs  -2.8  dB  for  non-parametric  baseline 

•  2-D  MSE-based  algorithm  was  .2  x  baseline  in  computations  and 
-.5  dB  better  than  baseline 
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BACKUP  #1  : 

ACTUAL  MCARM  RADAR  TEST  CONDITIONS 
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•  Performance  using  actual  MCARM  radar  data  with  injected  artificial  target  signals. 

•  Test  statistic  magnitude  vs  range  bin  index  for  acquisition  575,  flight  5  for  1  el  ch, 

4  az  channels,  32  pulses,  RBs  142-469  inclusive. 

•  Test  cases: 

•  Case  1:  fixed  windows  256  RBs  (142-269,341-468),  injected  SINR  -30dB  at  RB  291  and 

293;  plot  RBs  270-340. 

•  Case  2:  moving  window  11  RBs  (8  used);  otherwise  same  conditions  as  Case  1 

•  Case  3:  4  secondary  interfering  targets  of  same  doppler  added  in  RBs  238,  269,373,  400 

at  -10  dB  input  SINR 

•  Performances  not  sensitive  to  model  order. 

•  Measures:  max  target  peak  above  the  OdB  reference  level;  difference  between  target 
peak  value  and  highest  non-target  peak  value  (want  both  to  be  large) 

•  Prior  paper  found  AMF  poor  relative  to  PAMF. 

•  All  plots  normalized  to  common  mean  level  of  0  dB  to  set  reference  level 

•  LEFT  FIGURE:  2006  Study  ;  RIGHT  FIGURE:  2000  Study 
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BACKUP  #2  :  ACTUAL  DATA  TEST  CASE  3- 
MSE-BASED  M-C  AR  PARAMETRIC  ALGORITHM 
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J= 4  Array  Elements,  N=32  Pulses,  K= 256  Fixed  Range  Bin  Window, 
4  Interferers 


Sliding  Window,  Order=1 


Sliding  Window,  Order=2 


Range  Bin  Index 
Sliding  Window,  Order=3 


Range  Bin  Index 
Sliding  Window,  Order=5 


Range  Bin  Index 


Range  Bin  Index 
Sliding  Window,  Order=4 


Range  Bin  Index 
Sliding  Window,  Order=6 


Range  Bin  Index 


2006  STUDY  (orders  1  to  6) 


2000  STUDY  (order  4) 
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BACKUP  #3  :  ACTUAL  DATA  TEST  CASE  3- 
LSE-BASED  M-C  LP  PARAMETRIC  ALGORITHM 
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1=  4  Array  Elements,  N=32  Pulses,  K= 256  Fixed  Range  Bin  Window, 
4  Interferers 


Sliding  Window,  Order=1 


Range  Bin  Index 


Sliding  Window,  Order=3 


Range  Bin  Index 


Sliding  Window,  Order=5 


Range  Bin  Index 


Sliding  Window,  Order=2 


Sliding  Window,  Order=4 


Range  Bin  Index 
Sliding  Window,  Order=6 


2000  STUDY  (order  2) 


2006  STUDY  (orders  1  to  6) 
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